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ABSTRACT 


The potential drop in the polar cap region of a rotating, 
magnetized neutron star can be found assuming that the magnetic 
field is dipolar, with the field aligned (or anti-aligned) with 
the rotation axis. The curvature of the field lines is of critical 
importance. Charge flow is assumed to be along magnetic field 
lines. The electric field has a maximum at radius 1.5R* and the 
magnitude and functional form of the current can be determined. 
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I. Introduction 



The source of the accelerating potential in pulsars continues to 
be a significant problem in our understanding of these objects. It is 
clear that particles are accelerated but the source of the accelerating 
electric field is not yet understood. The most widely accepted picture 
of a pulsar is one in which particles are accelerated in the region of 
the magnetic polar cap, resulting (via some kind of plasma or MHD 
instability ) in a pencil beam of radio emission ( Radhakrishnan and 
Cooke, I 969 ) . Specific models along these lines have been proposed 
by a number of investigators (Pacini and Rees (1970)> Sturrock (I97I), 
Roberts and Sturrock (1972), Ruderman and Sutherland (I 975 )), 

In this paper we will investigate the form of the electric 
potential in the polar cap region of an aligned magnetic star. 


Clearly this cannot be a realistic model of a pulsar, but it can 
give us some insight into the physics of a real pulsar. We shall 
assume that the magnetic field is dipolar and the polar cap region 
is bounded by the first closed field line (i.e. the first line that 
does not connect to the surrounding nebula) . The field lines are 



of Sturrock (I 97 I) and Ruderman and Sutherland (1975)* R is 

max 


1 


the radius of the light cylinder, whereas the work of Roberts and 
Sturrock (1972) it is the "force-balance" radius. 

II. The Polar Cap Region 

In the region of the polar cap, the magnetic field is much larger 
than any expected electric field, so that it is reasonable to assume 
that charged particles will be tied to the magnetic field lines. In 
this case the current density is approximately proportional to the 
magnetic field strength 



If the stream flows in one direction only, e.g. outward, the charge 
density is given by p=j/v. We will consider the situtation where v a c 
and thus p*j/c. Clearly this is not true in the region immediately 
above the polar cap, but we will show that it is a good approximation 
when the distance above the cap is small compared with the overall 
dimensions of the system. We also note that in this situation the 
component: of j rotating with the star is negligible until the particles 
approach the light cylinder. In the Roberts and Sturrock model (1972), 
the polar-cap region extends only to the force-balance radius and hence 
the co-rotating component of j is always negligible. 

Poisson’s equation can now be written as 


2 

v $ = 



2 . 2 ) 


2 



At this point we introduce orthonormal dipole coordinates defined hy 


£ 


(cos e) 


sin 0 


(2.3) 


The curves of constant n are the magnetic field lines and curves of 
constant £ are the lines of constant magnetic potential. The remaining 
coordinate is the azimuthal angle <p . Since the system is azimuthally 
symmetric, we can now write Poisson’s equation as 


7** - l.+.3cgL-0 r^_ /I + 

4r 5 (cos 0)*s 


f— /— — 1 + (2.4) 

19£ \£ 9£ / n Zr)\ 3n /J 


where r and 0 are the usual spherical co-ordinates and are implicitly 
defined by £, n. 

We now make a major simplifying assumption, that the contribution 
of the divergence of the electric field due to the field component 
transverse to the magnetic field is small compared with that due to the 
field component parallel to the magnetic field. This is opposed to the 
assumption of Goldreich and Julian (1969), according to which the parallel 
component vanishes. We are implicitly assuming that the assumption of 
current outflow from the polar cap necessitates the development of a 
parallel electric field, and we are studying the effect of this component, 
whi<~h is associated with current outflow, independently of the transverse 
electric field component, which is associated with corotation. A complete 
solution of the electric field configuration of pulsar magnetospheres 
will require the simultaneous incorporation of both effects. However, 
since much work has been based on the Goldreich-Julian simplification, 
it seems not unreasonable to explore briefly the implications of the 
converse simplification. 
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We therefore now assume that the n derivatives can be neglected 


with respect to the £ derivatives; 

l+3cos z 0 9 /I 3*\ 47TJ * R * 


V' $ rs 


4r 5 (cos 0)«s 8£ 


/£ 

\£ 3£/ 


From equation (2.5), it immediately follows that 

•£ 


i a# 

£ n 


£ 4tt1 


£, 


,r- 

i 


(cos 0) 


d£ 


l+3cos 0 


The integration is along a field line n*const., and we can 


express d£ in terms of dr. Since 


d£ = - 


(cos 0)^dr 


and 


2r 


3 / 2 


sin 0 d0 


2r(cos 0) 


. sin 0 . 

dr) “ - - — i — dr + 


cos 0 


d0 = 0, 


we see that 


d£ = - 


l+3cos z 0 
4r z (cos 0) 3/z 


dr 


so that equation (2.6) becomes 


1 d$ 

£ dC 


R * 

c * 


f 


dr 


cos 


The equation of a field line is then used to write 


cos 0 * (1-ar) 


where a is a constant for each field line, given by 


( 2 . 5 ) 


( 2 . 6 ) 


therefore 


( 2 . 7 ) 


( 2 . 8 ) 


( 2 . 9 ) 


( 2 . 10 ) 
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( 2 . 11 ) 


a « sin 2 6./R 




The integration is now trivial: 


ff .[^ k [(1-o.R^ - (l-ar)^] + l |f |]t 


( 2 . 12 ) 


We may integrate equation (2.10) again, using the condition that we are 
integrating along a field line, and so obtain 




ca 


[fc (b - 7 ) - ~ (r - r;) + 


2 _ 2 

, i U d± 

2 df 


(2.13) 


+ ♦(«:*) 


At this point we note that the dependence of <5> on n (or angle 0*) is 
completely determined by specifying $ and d<J>/d£ at £ * £* as functions of n. 
We also note that the value of £* is itself a function of n (i.e. of 
the particular field line chosen). However, one may find a surface 
f a = £ which is very close to the surface of the star in the polar-cap 
region. Since a neutron star is a nearly perfect conductor, the force 
on a charged particle rotating with the star must be zero. This means 
that, on the surface of the star. 


E a = -(:: X r)/c X B, 


(2.14) 


which in turn implies that on the stellar surface 


2c 


(cos' , 0. - cos ? 0 


(2.15) 


5 


We have used our freedom to change t by a constant to set on the 
closed field lines that bound the polar-cap regions. 

If we now require that the charge flow be space-charge limited 
(this condition has been used in a number of models, e.g. Sturrock (1971), 
we require that the electric field parallel to the magnetic field lines 
be zero at the stellar surface. This implies 


d$ 

d£ 


( 2 . 16 ) 


Finally, if we can determine either $> or d$>/d£ at an outer boundary 
£ ■ £ , we can determine j # . Specifically, we note that d4>/d£ i 0 

for any finite distance above the stellar surface. 

To see more clearly what these results mean, we can look at equation 
(2.12) in the limit Q* « 1: 


4> s 


8 tt J* 

c 




(2.17) 


Thus we see that in this limit 4> appears to behave in the same way as 
the free-space electric field produced by a monopole and a dipole. We 
note also that the maximum electric field occurs at r - 3 /2 R*. 


III. Two-Dimensional Approximation 

In this section we shall attempt to improve the model by considering 

the two-dimensional case. The potential <t> can be expanded in a power 

series in n* The symmetry of the problem guarantees that the series has 

only even powers of n. Since n * 1/R ^ , n may be treated as 

small and we diall consider the first two terms only in the expansion. 

The boundary condition on the field lines n * n is t ■ 0, which 

max 
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determines the expansion: 


4> ( ' » n ) * $(C, 


7 n 2 ; 

max 


(3.1) 


On substituting this expression into equation 2.4, we obtain 


72 , . !+j£2 sfe L (l . IttiSi (l+3cof* 28) 

4r5(cos0)^ U 9K I r 3 r 2 ' ***> 

r 'max 


(3.2) 




Wr may now choose the field line n=0 on which to solve equation (3.2) 
On this field line 0=0 and r=l/£ so that 


_ (l _ 

K H ;if ./ “ 


4 X 3 4 « . R 3 c 3 

n ^ * 


(3.3) 


A particular solution to the inhomogeneous equation is 


■t (C,0) * R] t\ 2 


(3.4) 


The homogeneous equation can be solved by making the substitution 


u = f 2 and 4-a ,0) = u~B(u,0) 


(3.3) 


The function B(u) must then satisfy the equation 


^ + 1 ^-r 4 - + ^-l b-o 

du u du L u' n ■ J 


(3.6) 
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This is the modified Bessel equation of order 2 and parameter 4/n . 

We now have 


♦ (£, n) 



+ c 2 £K 2 





+ - 

c 




2 

max 


1 (- n4- ). 

J \ max f 


(3.7) 


We require that $ 
this condition. 


0 as f, + 0 (r * 

As z K), K 2 (z) ' — 
Z 2 


*). The I 2 term automatically satisfies 
so that we have 


♦( 0 , 0 ) 


C 2 n h R * 

8 n max + c 


max 


(3.8) 


which gives 


C 7 


At 5 



we want 


ne*,o) 





'*B* R * 


2c 


n 


?. 

max 


(3.9) 


(3.10) 


Since >>1, we can use the asymptotic expansions of I 2 and K 2 

to determine the constant cj. 


12 


z 


e 

(2TTZ) 13 


e Z with z 


wT 

n max 


Hence 

c 1 “ 



n '' 

max 




max 





( 2 tt z ) * e 


-z 


(3.11) 


(3.12) 
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If we now require the flow to be space-charge United, we can determine 
We write d$/d€ in terns of z (as defined in equation 3.11): 


<»♦<£. Q) . 8 Ii £[<1,2 l2 (1) 1 + 

dC zn 2 dz J.16 “** 2 1 J 


max 


p, 4r[w «»<■>] 

n max 


Using 


^ [z 2 I 2 (z)] - z 2 Ij (z) 


we find 


d£a,0) 

dC 


jjz 2 I,(.) 


Tz * 2 K i (l) 


Thus 


-2z 

Cj ■ +*e c 2 


- 8tt 2 J* R 3 e -2^ 


Combining (3.16) with (3.12) then gives us the following 
on j*: 

■ (f) 4 

Since z>>t, we can approximate j* by 


J* 


flBJLn 2 


\2n / 


* max 


(2ttz) 


max 



(3.13) 


(3.14) 


. (3.15) 


(3.16) 

condition 


(3.17) 


(3.18) 


which is the value estimated by Sturrock (1971). 



We would also like to know how varies across the surface of the 
star. To determine this, we return to equation (2.4) and write it as an 
expansion in powers of n. From the lefinition of £ and n we can write 


£ 4 r 4 + r n 2 - 1 


(3.19) 


We then solve for r implicitly, giving us 


r ■ 


' rW 




(3.20) 


This becomes (to ord n ” n 2 ) 


H* • «) - ° ( r, %4 <3 - 2i > 


This expansion is valid as long as n 2 < 4£ (i.e, as long as 6 £ 85°). 
We also must expand the potential to order n 2 giving us 


m,n> - *<£,0) [l - + a(On 4 ] + 0 (n 6 ) 

*■ n max ■* 

and we expand the current density to order n 2 giving 
J*(n) - j*(0) (1 + dj n 2 ) + 0(n 4 ) 

Equation (2.4) then becomes 

(*♦ lr) \h (1 $?“’ 0> )(> - ^-) 

1 ’ max' 


; 2 (- h +(vrr + «•«>) -> 2 ] 

£- *• max ' max ' 

4uj.(0) i t 2\ 

- __ R 3 t 3|, + li.) (1 * dj n 2 ) . 


(3.22) 


(3.23) 


(3.24) 
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V*.e zeroth order equation is equation (3.3) while the second-order 
equation is 



(3.25) 


Using the zeroth order equation to simplify equation (3.25), we obtain 

4"j*(0> R l 


- «♦«,<» + 


n *♦ 

max 


n 2 

max 


+ (t^V + 16a ^>) *<e.°) = 

' max / 


(3.26) 


- ^j*(0) 


*2 d , 


Since the variation of j with respect to r has already been included to 
all orders ( j -1/ r 3 ) we know that dj is independent of £• Equation 
(3.26) is satisfied for any form of 4> if 


a(0 


4n 2 
max 




and 


(3.27) 


d . — 
max 


(3.28) 


Thus 


j*(n) = - 


£2B 

2 




(3.29) 


and we have also determined the corrections to $ of order n 14 . 
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The fact that has the same n dependence (to second order) as 
4> suggests that we may be able to find other solutions of Poisson's 
equation for which this Is true. We therefore write equation (2.4) 
in the small angle approximation under this assumption. We write $ 
and as follows: 


*(S»n) - *(5)H(n) 

j*(n) * j*H(n) 


(3.30) 


In the small angle approximation Poisson's equation now becomes 

CH(n) - H(n)tf>*(£) + i|>(S)H"n + -*U)H 1 (n) 

n 

(3.31) 

- — R jH(n) . 

Dividing by we obtain two ordinary differential equations 


H”(n) + ^ H^n) + A 2 H(n) - 0 (3.32) 

4irj* 

« *"(5) - ^(e) - A 2 I KC) - “ Rl (3.33) 

c * 

The solution of equation (3.32) is 

H(n) * J 0 (An) (3.34) 


while the solution of equation (3.33) is 


*tt) m c.?I,0£) + c,CK ,(*£) + 




A z c 


Rl 


(3.35) 


The boundary condition on the closed field lines implies that 

-V x, w> - 0 «•»> 
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To complete the solution, the constants c 1 and c 2 must be determined 
by the boundary conditions at the stellar surface and at Infinity. This 
requires that be a sum of the form: 


<I>U) 


+ c u SK 2 (a/<) +I7r i R 
1 L 


(3.37) 


where the are determined by the zeros of equation (3.36) and 


j*(n) 


^ ^i J o (A i n) 


(3.38) 
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IV. Conclusion 


We have determined the analytic behavior of I’ by several different 
methods. The firm conclusion of these various approaches is that the 
current density is of the order ~ at the surface. The one 
dimensional approximation indicates that the electric field has a maximum 
on the order of one star radius above the surface. We have shown that 
it is possible to find a two-dimensional solution which satisfies the 
boundary condition on the star, goes to zero at large distances and 
results in space-charge limited current flow. 

We are currently preparing a numerical analysis of the two-dimensional 

solution in order to plot the surfaces of constant electric potential 

and determine the location of the maximum electric field for various 

values of n , B., and Q. 
max * 
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